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ABSTRACT 


This is the final report for NASA Contract NAS 2-8664, "An Analysis of Water in 
Galactic Infrared Sources Using the NASA Lear Airborne Observatory. " 

We have completed a very successfid program using the Grumman interferometer 
system on the NASA Lear Jet Airborne Observatory. Following a description of the instru- 
mentation and data reduction procedures, a review is presented of the objects observed 
(standard stars, M stars, a nebula, planets, and the moon). The observing parameters are 
listed for each flight date. The spectra obtained from these data flights are presented, 
grouped by class of object. A brief discussion is presented of how these spectra can be 
utilized. 

A number of presentations and publications have already resulted from this contract 
and are listed in this report. 
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INTRODUCTION 


This is the final report for NASA Contract NAS 2-8664, "Analysis of Water in 
Galactic fiifrared Sources Using the NASA Lear Airborne Observatory, " from 2 January 1975 
to 31 March 1978. We describe the instrumentation used, the observational data obtained, 
and our analyses to date as well as the potential for future work. 

Observational data were obtained on every scheduled flight, and met or exceeded the 
contract requirements. In addition, an order of magnitude better signal -to-noise ratio (S/N), 
was obtained than originally anticipated allowing the full resolution capability of the system to 
be realized. We report here in detail on those spectra with high S/N with all other spectra 
being available on request in a tape format. We have analyzed some of the surprisingly large 
quantity of excellent data, presented this material orally at meetings in the scientific com- 
munity and published the detailed results in the open literature. (See Publications and 
Presentations section.) 

The objects observed are standard stars, M stars, a nebula, planets, and the moon. 

The M stars have large amoimts of water in their atmosphere and can provide information 
about the evolution of water in the galaxj'-. The emission lines in the nebula (M 42 in Orion) 
are useful in determining the reddening and dust properties in the ionized region. The 
planetary observations reveal information on the physical properties of the atmosphere and 
surface, and the lunar observations observed through the earth's atmosphere provide infor- 
mation on the chemical composition of the stratosphei’e. The standard stars provide calibra- 
tion checks on the above objects. 

INST RlIM E NTAT ION 

In the field of infrared astronomy the quality of the instrumentation used has always 
been a strong factor in determining the success of an experiment. This requirement is 
especially strong for airborne infrared astro lomy. We report here on a general description 
of the Grumman Michelson interferometer instrumentation, and in particular, the design of 
the optics, and the use of advanced data reduction procedures. 

General Instrumentation 

The Michelson interferometer remains the ideal instrument to be used for the observa- 
tion of water vapor in M stars and the interstellar medium. Its large throughput (no slits) and 
multiplexing ability (all wavelengths measured simultaneously) allows miiximum efficiency 
when observing faint objects (K mag 0) at good wavelength resolution (2 cm~^) over a large 
bandpass (1.5-3. 5 um). An entire spectrum is obtained in one 2-second scan with the S/N 
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inipTOved by means of co-adding hundreds of scans. Guiding problems are thereby reduced 
over those required by a scanning spectrometer. 

The Lear Alrbonie Observatoiy, flying at 13.7 km (above the tropopause), reduces 
the overhead precipitable water vapor to less than 10 ^-m compared to several mm .ver a 
ground based observatory. The availability of frequent observing flights with the Lear Jet 
allows us to observe how the properties of the M stars vary with phase. 

The configuration we use consists of a NASA heliostat and a vibration-isolated 30-cm, 
f/23 telescope hard-mounted to the interferometer assembly. Figure 1 shows a photograph 
of the system mounted in the NASA Lear Jet. On the left is the heliostat mirror which com- 
pensates for the airplane motions. The starlight comes in through the infrared transmitting 
cabin window (not shown) and is reflected off the heliostat into the telescope shown in the 
center and extending to the right. Further to the right and lower in the picture is the mount- 
ing platform containing foreoptics, dewars, and the interferometer. Figure 2 shows this 
same optical system schematically with more detail on the light path through the foreoptics 
and interferometer. The foreoptics is of a symmetric, dual beam type; onl 3 '^ one side of the 
light path is shown traced to the detectors. The telescope is a folded Dall-Kirkham type 
designed by us to give a field of observation of 6,3 arc-minutes. This wide field is a neces- 
sary requirement to observe the extended object in our program, M42 (the Orion Nebula). 

The Grumman interferometer system operates in a double beam, double detector configura- 
tion to give automatic cancellation of the emitted light due to the sky, telescope, and inter- 
ferometer. The detectors are nitrogen-cooled 1/2 mm LiSb cells with cold filter limited 
response out to 3.5 ^m, where the GE 105 window effectively stops transmitting. The 
interferometer is set to run at any resolution down to 2 cm which is the maximum possible 
with an f/23 beam. The telescope interferometer system floats on a foam pad resting on a 
plate bolted to the baggage compartment floor. The system has been developed to the point 
where both mechanical and electrical interference from the aircraft have been overcome. 

Our guiding telescope was initially capable of locating objects as faint as V = 10. This 
guider has been upgraded by the addition of an image intensifier. The visual magnitude limit, 
since we are looking through the main telescope, is now V = 12, which w'as sufficient for all 
our observations except for a few of the M stars at minimum phase. For all these observa- 
tions offset guiding was used. 

Data from the interferometer were recorded on a instrumentation grade analog tape 
recorder and digitized on the ground after every flight. The digitization and co-addlng was 
done on the Grumman Research Department's Portable Data Acquisition System which also 
allowed data to be sent by phone link back to the Depai'tment's IGK word computer. A 
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complete full resolution spectrum was computed and sent back to us at Ames for initial 
evaluation. We were then able to use this information in planning the observations for the 
next flight. For example, the interferometer mirror travel (which determines spectral 
resolution), and total observing time (which determines S/N) occasionally had to be changed 
based on previous flight data. One of the strong limitations of IR interferometry in the past 
has been the inability to look at the spectrum in "real time. " Usually, the astronomer only 
sees a spectrum after he has left the observatory and finished his series of measurements. 
Our data system and the necessary computer programs are developed and have been used 
successfully for our four flight series. 

Optical E)esign 

An important source of progress in our optics laboratory has been use of a computer 
optical design capability, allowing us to design foreoptics to mate the interferometer to var- 
ious telescopes in a way that is optimum for the type of object to be observed. Figure 3 
shows a typical output from this program in which it was desired to fit the interferometer to 
an f/23 telescope with a 6.3 arc-minute field of view. In this case, available room was also 
tightly constrained so that the beam had to be folded several times. Figure 4 shows a picture 
of this foreoptics system in our laboratory in preparation for testing with our interferometer 
and detectors. In this system, which is the one eventually used in the Lear, we have also 
included remotely adjustable alignment mirrors and focal plane guiding optics for ease of 
operation in an aircraft. We have used this system to make low resolution observations 
from the ground at our own laboratory, using a 30 cm highly folded telescope we designed and 
built at Grumman. Figure 5 shows an example of our results on a humid night. In addition 
to the atmospheric absorption by CO 2 , one can see the very strong water absorption features 
which confuse groundbased astronomical observations of water. 

Data Reduction Procedures 

After every flight the tape recorder was taken from the plane to the computing facility 
brought to the Ames Research Center. The tapes were first played through, and all verbal 
comments were transcribed along with a relative time code which was recorded on the tape 
prior to each flight. This code allows reliable repeat access to any region of the data and 
was found to be much more accurate than the conventional tape recorder foot counter. 

The tape recorder data channel was connected to a Hewlett Packard (HP) 5465A A/D 
converter and the white light (start) and HeNe laser (sampling) channels were connected to a 
GAC -designed interface electronics unit which generates synchronized pulses and sends them 
to the A/D unit. Synchronization with the desired data is thus assured independent of any 
tape recorder or interfeometer speed fluctuations. 
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A Texas Instrument Silent 700 Terminal with dual cassette reader was connected di- 
rectly to the HP 21MX computer (a minicomputer) and the 21MX programs BABL8 and 
LEWDAL were loaded into the computer. These programs (stored on tape cassettes) were 
previously generated by the Grumman Research Department and allow the HP 21MX to 
interact \vith the larger time-sharing HP3000 computer located in the Research Department 
Computing Center in Bethpage, New York. From discussions w^th HP, this represents the 
first successful attempt of a direct link between an HP3000 and an HP 21MX. 

A junction switch provides a TI terminal-HP3000 or HP 21MX-HP3000 connection 
through the phone lines. This transcontinental computer linkup has worked successfully 
and allows each night's data to be examined the next day. In the Jan/Feb 1976 series of 
flights, after the first 5 flights had demonstrated the good quality of the data we were obtain- 
ing, we decided to increase the spectral resolution from 100 to 13.3 cm"^. Without the 
ability to see the finished spectra, we would have been more conservative in changing our 
resolution. 

A BASIC A/D program, which determines the number of interferogram scans to be 
co-added and the number of points to be digitized in each scan, is transmitted from the 
HP3000 in Bethpage to the HP 21MX at Ames. The data tape is started at the desired scan 
and the digitization proceeds. Scans with the star in one aperture are co-added and the re- 
sult subtracted from the co-added scans with the star in the other aperture. In this way the 
residual background from the sl<y, airplane, telescope, and interferometer is considerably 
reduced. The individual scans are visually monitored on an oscilloscope and any that are 
considered irregular can be knocked out of the co-adding process through a switch on the 
HP 21MX. After co-addmg, the resultant digital interferogram is displayed on the HP 5460A 
display unit. If the interferogram has sufficient S/N it is transmitted in a packed format over 
the phone lines to the HP3000. The transmission takes about 4 minutes for 7100 points at 
1200 BAUD. If for any reason a "glitch" develops in the transmission due to noise on the 
long distance lines, the data can be retrieved by the computer and retransmitted without 
having to redigitize the data. 

All operations on the data are now done by the HP 3000 through the TI terminal. The 
data in the HP 3000 is unpacked (symbols converted to digital numbers) and the region near 
the interferogram white light pobit is transmitted back and plotted on the terminal at Ames. 

To maximize observational efficiency, the interferogram is generated single-sided, i.e. , 
the white light point is to one side of the data. This increases the S/N by for a given 
observing time, but the interferogram needs to be phase corrected before the Fourier 
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transform can be taken. The phase correction technique used is that if Forman et al (1966) 
and the correction program PHASE is listed iri the Appendix. 

This program symmetrizes the interferogram about the white light point; hence, a 
cosine transform can be done on the data. A fast cosine transform program HFFT (see 
Appendix) using the Cooley-Tukey method has been developed at GAC which is extremely 
efficient in computer core storage, i.e. , the transform of S192 points only needs 8200 stor- 
age locations. This program also smooths the spectrum with a Hamniing function (Blackman 
and Tukey, 1958) and interpolates for four points between each pair of output points for ease 
in line identification. To our loiowledge, no fast cosine transform program has previously 
been published which is even within a factor of two of the storage efficiency of this program. 

The spectral resolution is given by 

/ 7901.14 \ , 

Av = 1.815 * ( 1 cm” 

where N is the number of data points from the white light point, 7901.14 cm ^ is the Nyquist 
frequency when the He-Nc laser line at 0.63282 urn is used as the interferogram digitizing 
clock (7901.14 cm”^ = 1/ [2 x .63282 x 10 cm] ) and 1.815 is the rai'o of the full-width 
at half maximum (F\VHM) of the Hamniing function to the frequency separation of the output 
spectral points (7901. 14/N). Tliis frequency separ^xion also equals the first zero of the sine 
function (sin(x)/x) which is the inherent apodization of the spectrum with no smoothing. The 
familiar "ringing" phenomenon is suppressed by convolving the spectiaim with the Hamming 
function (Bell, 1972). 

Selected regions of the complete spectrum, i.e., CO 2 bands at 3700 cm”^, are 
transmitted back to the TI terminal and plotted, and the S/N can be examined to see if it is 
sufficient. Since all the reduction procedures just described can be done after each flight, 
necessary changes to the equipment or experimental procedures can be made for subsequent 
flights. This real time look at the finished sp>ectra has proven invalualile to the success of 
our program. Upon completion of the flight series and return home, all s;x'ctra arc 
plotted on a Tektronix 4014 Graphics teiminal and a hard copy is made. 

OBSERVATIONS 

An observational data summaiy is presented in Tables lA and IB for tlic 36 flights 
conducted in 1976 and 1977. Column 1 is the flight number, 2 is the flight date in U.T., 3 
lists the objects observed, 4, 5 and 6 are the interferometer settings (ix;lativc), 7 is tlie 
data filter in Hz, 8 is the amplifier gain in db with the first value referring to the first ob- 
ject, etc., 9 identifies spectra shown in the report, and 10 lists particular comments 



explained at the bottoni of Table IB. All observations used a 6.3 arc-min aperture ex- 
cept for obser\'ations of the moon after Flight 17. 

The primary goal of our progi’am is its scientific return. The spectra shown are the 
first step in obtaining this return, and are in compliance with the requirements of the con- 
tract work statement. A brief discussion of some of the results follows. A list of spectra is 
given in Table 2 arranged by class of object-standard (non -variable) stars, M stars, nebula, 
planets, and stratospheric observations. Column 1 is the object observed, 2 is the flight 
number, 3 is the date in U.T. , 4 is the number of scans, 5 is the spectral resolution in cm ^ 
usmg Hamming apodization, 6 is the file name used for identificati purposes by the com- 
puter, and 7 is the figure number m this report. 

For observations of the moon (diameter ~ 30 arc-minutes), tw'o separate apertures 
(labeled L and R) were used and these are so labeled in Table 2 and in the figure. 

The w'avenumber range shown m all spectra is 2800-6000 cm ^ with the limits being 
determined by the optical filter m front of the InSIj detectors. 

The captions on all figures identify the object, the date of observ'ation (U.T.), and 
the resolution in cm ^ defined as the full-wddth at half maximum of the Hamming apodization 
function. For some of ihe spectra the resolution has been degraded from the maximum 
measured values to improve the S/N value. 

Per the w'ork statement, these spectra are shown in relative units uncorrected for 
frequency shift, mstrument response, and residual telluric absorption. As such their pri- 
mary ability is for spectroscopic identification of molecular features. Hown'-er, w'e have 
reduced a number of selected spectra further for purposes of publication. An example of the 

«og 2 

technique used to reduce these spectra to absolute fluxes ui Janskys (1 = 10 “ W/m /Hz) 

and true (laboratory) frequencies is discussed later in this report along with several 
samples of sjjectra so corrected. 

DISCUSSION OF SPECTRA 

The results of the observations presented ui the previous section arc now discussed 
accordmg to class of object -standard siars, M stars, nebula, planets, and stratospheric 
observations. A brief discussion is given on important features seen in the spectra and on 
their significance. This is followed by examples of reduction to absolute flux. 

Standard Stars 

The standard (non-variable) stellar sijectra are shown in Figures 6 - 23 for a Ori, a Tau, 
a Boo, y. Cep, a Lyr and a CMa. The only stellar features easily seen arc the first 
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TABLE lA 1976 OBSERVATIONAL DATA 


Relative Literferometer Settings 
(4) (7) (8) 


(1) 

(2) 

(3) 

Sweep 

(5) 

(G) 

Data 

Ampl 

(9) 

(10) 

No. 

Date (UT) 

Object 

Length 

Velocity 

Offset 

Filter (Hz) 

Gain (db) 

Spectra 

Comments 

1 

27 Jan 

Moon, Venus 

0.010 

100 slow 

840.5 

50-100 

15,40 


A 

2 

28 Jan 

O' C Ma, or Ori 

0.010 

100 s 

840.5 

50-100 

35,35 


A 

3 

29 Jan 

RX Boo, a Tau 

0.010 

100 s 

840.5 

50-100 

30,30 



4 

30 Jan 

MCep 

0.010 

100 s 

840.5 

50-100 

30 



5 

3 Feb 

»C Ma, orOri 

0.010 

100 s 

832 

50-160 

30,25 



6 

5 Feb 

OfOri 

0.040 

195 fast 

750 

100-315 

35 

<sr Ori 


7 

6 Feb 

RX Boo, (yC Ma, 
or Tau 

0.040 

195 f 

750 

100-315 

45,45,40 

RX Boo, or Tau 


8 

10 Feb 

01 C Ma, Moon 

0.040 

195 f 

750 

100-315 

45,10 

Moon 


9 

11 Feb 

Mars, Moon 

0.040 

195 f 

750 

100-315 

50,10 

Mars, Moon 


10 

12 Feb 

Jupiter, Mira 

0.040 

195 f 

750 

100-315 

40,40 

Mira 


11 

24 Aug 

a Boo 

0.020 

128.5 f 

420 

40-160 

35 


A 

12 

25 Aug 

a Boo, RX Boo 

0.270 

370 f 

120 

250-1000 

35,35 



13 

26 Aug 

a Boo, RX Boo 

0,270 

370 f 

120 

250-1000 

40,40 



14 

27 Aug 

a Boo, RX Boo 

0.270 

370 f 

120 

250-1000 

45,45 

O' Boo, RX Boo 


15 

31 Aug 

a Lyr, M42 

0.020 

128.5 f 

420 

40-160 

45,45 



16 

1 Sept 

U Ori, a Ori 

0.270 

370 f 

120 

250-1000 

55,45 

a Ori 


17 

3 Sept 

Moon 

0.270 

370 f 

120 

250-1000 

15,45 

Moon 


18 

8 Sept 

a Lyr, M42 

0.047 

185 f 

391 

80-315 

45,45 


B 

19 

9 Sept 

or Lyr, M42 

0.055 

185 f 

384 

80-315 

45, .5 

O’ Lyr, M42 

B 

20 

10 Sept 

Of Boo, Moon 

0.270 

370 f 

123 

250-1000 

45, 25 

Moon 

B,C 
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TABLE IB 1977 OBSERVATIONAL DAIA (Cont) 
Relative Interferometer Settings 


(4) (7) (8) 


(1) 

No. 

(2> 

Date (UT) 

(3) 

Object 

Sweep 

Length 

(5) 

Velocity 

(6) 

Offset 

Data 

Filter (Hz) 

Ampl 
Gain (db) 

(9) 

Spectra 

(10) 

Comments 

21 

8 Feb 

CMa, y Ori 

055 

185 f 

384 

80-315 

50,40 

a CMa, a Ori 


22 

9 Feb 

ot CMa, 

055 

185 f 

384 

80-315 

45 

a CMa 




Moon 

270 

370 f 

123 

250-1000 

20 

Moon 

C 

23 

10 Feb 

Cep, Mira 

055 

185 f 

384 

80-315 

45,40 

p. Cep, Mira 


24 

11 Feb 

a CMa 

055 

185 f 

384 

80-315 

45 

a CMa 


25 

17 Feb 

RLeo, Venus 

055 

185 f 

384 

80-315 

45,30 

RLeo, Venus 


26 

18 Feb 

RLeo 

055 

185 f 

384 

80-315 

45 

RLeo 


27 

12 Oct 

UOri, or Ori 

055 

185 f 

384 

80-315 

40,35 

ot Ori 

D 

28 

13 Oct 

UOri, a Ori 

055 

185 f 

384 

80-315 

50,40 

a Ori 

D 

29 

14 Oct 

UHer, Mira 

055 

185 f 

384 

80-315 

45,45 

Mira 

D 

30 

18 Oct 

RLeo 

020 

228,5 f 

429.5 

100-315 

40 

RLeo 

D 



u Cep 

055 

185 f 

384 

125-315 

40 

p- Cep 

D 

31 

19 Oct 

UHer 

020 

228.5 f 

429.5 

100-315 

50 

UHer 

D 

32 

20 Oct 

Moon 

270 

370 f 

123 

250-1000 

30 

Moon 

C,D 



RAql 

020 

228.5 f 

429.5 

100-315 

40 

R Aql 

D 

33 

21 Oct 

UHer 

020 

228.5 f 

429,5 

100-315 

45 


D 



Mira 

055 

185 f 

384 

125-315 

45 


D 

34 

26 Oct 

a Ori 

020 

228,5 f 

429.5 

100-315 

45 

ot Ori 

D 

35 

27 Oct 

a Ori, M 42 

020 

228,5 f 

430 

80-250 

40,45 

ot Ori, M42 

D 

36 

28 Oct 

Mars 

055 

185 f 

385 

125-315 

50 

Mars 

D 



ot CMa 

020 

228,5 f 

430 

80-250 

45 

ot CMa 

D 


Comments: 

A. Dewai D 11 (rear) only 

B. Block Jig filter with long wavelength cutoff of 3.15 um used. All other flights used filter with cutoff of 3.55 um. 

C. One arc-minute limiting aperture placed in front of Dewar D 10 (front) for moon observation to increase spectral 
Resolution. Dewar D 11 (rear) was disconnected. 

D. Observations were conducted using a quartz window (S/N 016) which contained water of hydration bands. 



TABLE 2 LIST OF SPECTRA - STANDARD STARS 


(1) 

object 

(2) 

Fit 

(3) 

Date (UT) 

(4) 

No, Scans 

(5) _i 

Resol. (cm ) 

(6) 

File Ntunc 

(7) 

Fig. No. 

a Ori 

G 

5 Feb 197G 

G4 

12.3 

AORIGSP 

6 a, b 

Of Ori 

IG 

1 Sep 197G 

170 

1.8 

AORlloSP 

7a, b, c, d 

01 Ori 

IG 

1 Sep 197G 

100 

2.1 

AORIIGSI 

8a, b, c, d 

Of Ori 

21 

8 Feb 1977 

250 

9.9 

AOR121SP 

9 a, b 

a Ori 

27 

12 Oct 1977 

200 

9.9 

AORI27S1 

lOa, b 

ot Ori 

28 

12 Oct 1977 

140 

9.9 

AORI28SP 

11a, b 

O' Ori 

24 

2G Oct 1977 

100 

27.7 

\ORI24SP 

12 

a Ori 

25 

27 Oct 1977 

150 

27.7 

AORI25SP 

13 

Of Tau 

7 

G Feb 197G 

140 

12.2 

ATAU7SP 

14a, b 

o Boo 

14 

27 Auk lf>7G 

170 

1.8 

AB0014S1 

15a, b, c, d 

a Boo 

14 

27 Aug 197G 

100 

2.1 

AB0014SP 

IGa, b, c, d 

P. Cep 

22 

10 Feb 1977 

200 

29.9 

MCEP22S2 

17 

p Cep 

20 

18 Oct 1977 

100 

29.9 

MCEP20SP 

18 

o Lyr 

19 

9 Sep 197G 

500 

5(). 0 

ALY19BSP 

19 

ot CMa 

21 

8 Feb 1977 

180 

29.9 

ACM21S1 

20 

a CMa 

22 

9 Feb 1977 

200 

29.9 

ACM22S2 

21 

a CMa 

24 

U Feb 1977 

120 

29,9 

ACM24S1 

22 

o CMa 

2G 

18 Oct 1977 

100 

27.7 

ACM2GSP 

22 




M STARS 



RX Boo 

7 

G Feb 197G 

IGO 

29.9 

RB007SP 

24 

RX Boo 

14 

27 Aug 197G 

200 

29.9 

RB0014S2 

25 

Mira 

10 

12 Feb 197G 

290 

12.2 

MIRAIOSP 

26a, b 

Mi"a 

22 

10 Feb 1977 

200 

9.9 

M1RA22SP 

27a, b 

Mira 

22 

10 Feb 1977 

150 

12.2 

MIRA22S1 

28a, b 

Mira 

29 

14 Oct 1977 

100 

29.9 

M1RA29SP 

29 

R Leo 

25 

17 Feb 1977 

100 

9.9 

RLE025SP 

30a, b 

R Leo 

2G 

18 Feb 1977 

280 

9.9 

RLE02GS1 

21a, b 

R Leo 

2G 

18 Feb 1977 

200 

9.9 

RLE02GS2 

22a, b 

R Leo 

2G 

18 Feb 1977 

200 

12.2 

RLE02GS2 

33a, b 
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TABLE 2 (Cont) 


(1) 

Object 

(2) 

Fit 

(3) 

Date (UT) 

(4) (5) 

No. Scans Resol (cm ) 

(6) 

File Name 

(7) 

Fig. No. 




M STARS (Cont) 



R Leo 

30 

18 Oct 1977 

350 

37.7 

RLEO30SP 

34 

U Her 

31 

19 Oct 1977 

150 

37.5 

UHER31S1 

35 

R Aql 

32 

20 Oct 1977 

350 

37.3 

RAQL32SP 

36 





NEBULA 



M 42 

19 

9 Sep 1976 

500 

56.0 

M4219ASP 

37 

M 42 

IS 

9 Sep 1976 

500 

56.0 

M4219AB 

38 

M 42 

35 

27 Oct 1977 

300 

56.0 

M4235S1 

39 





PLANETS 



Venus 

25 

17 Feb 1977 

250 

9.9 

VEN25SP 

40a, b 

Mars 

9 

11 Feb 1976 

288 

13.3 

MARS9SP 

41a, b 

Mars 

36 

28 Oct 1977 

100 

29.9 

MARS36S1 

42 



STRATOSPHERIC OBSERVATIONS 



Moon (L) 

8 

10 Feb 1976 

30 

13.3 

MOON8LSP 

43a, b 

Moon (R) 

8 

10 Feb 1976 

15 

13.3 

MOON8RSP 

44a, b 

Moon (L) 

9 

11 Feb 1976 

15 

13.3 

MOON9I.SP 

45a, b 

Moon (R) 

9 

11 Feb 1976 

15 

13.3 

MOON9RSP 

46a, b 

Moon (L) 

17 

3 Sep 1976 

40 

1.8 

MOON17L 

47a, c, e, g 

Moon (R) 

17 

3 Sep 1976 

50 

l.S 

MOON17R 

47b, d, f, h 

Moon (L) 

20 

10 Sep 1976 

15 

2.1 

MOON20L1 

48a, c, e, g 

Moon (R) 

20 

10 Sep 1976 

15 

2.1 

MOON20R1 

48b, d, f, h 

Moon (L) 

22 

9 Feb 1977 

18 

2.1 

MOON22L1 

49a, c, e, g 

Moon (R) 

22 

9 Feb 1977 

26 

2.1 

MOON22R1 

49b, d, f, h 

Moon (R) 

32 

20 Oct 1977 

100 

1.8 

MOON32R 

50a, b, c, d 
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overtone bands of and (4000-43G0 cm~^) for the four late type stars. Other 

stellar features are much wealcer. The spike in the a Lyr spectrum (Fig. 19) at 3250 cm ^ 
is a result of uncancelled background. The a CMa spectra (Fig. 20-23) are noisy but provide 
a check on our calibration techniques. 

hi Stars 

The hi (vai'iable) stars are shown in Figures 24-30 for RX Boo, Mira, R Leo, U Her 

and R Aql. Strong, broad steam bands centered at 3700 and 5300 cm~^, as well as the first 
12 16 I 3 l6 1 

overtone bands of C 0 and '^O (4000-4360 cm ) are seen in all stars. 

Using the best available stellar model atmospheres and a line parameter atlas of high 
temperature steam, information about the temperature, water abundance, and surface gravity 
of the stars can be extracted. 

Nebula 

The bright nebula hI42 in Orion is shown in mlative units in Figures 37 and 39 and in 
Janskj’s in Fig. 38. The absolute curve for hl42 wns obtained by dividing the relative hl42 
curve (Fig. 37) by a Lyr taken the same night and at the same air mass (Fig, 19) and then 
multiplying by the Imown flux of a Lyr. The K magnitude at 4545 cm ^(2.2 um) is +1.7 
with a 6.3 arc-min FOV centered on Ori. The M42 data of Fit. 35 (Fig. 39) is noisy and 
was not reduced to absolute units as no suitable calibration object was observed during that 
flight scries (Oct. 1977). The strong hydrogen emission lines of Paschen a (5333 cm~^), 
Brackett p (3809 cm ^) and possibly Brackett y (4618 cm ^) ai’c seen and identified in the 
figures. 

Planets 

The two planets Venus and Mars are showai in Figures 40-42. Strong CO, bands at 
4850, 4980, and 5110 cm ^ and the w’uig of the very strong band at 3750 cm~^is seen on 
both planets, v/ith Venus showing additional bands at 4600 and 5300 cm”^. 

Modcliiig of the Martian atmosphere permits the extraction of Martian soil reflectance 
in the CO., bands centered at 3657 cm~^. We have concluded that the lUartian soil surface 
near-infrared optical propei'ties are consistent with a soil composition similar to montmoril- 
lonite or limonite, mixed with a basalt Egan et al, 1978). Also, it w-as foimd that the non- 
dust storm aerosol conditions existing m the Martian atmosphere can influence COg line 
shapes and, hence, the average transmission. Thus, planetary spectroscopic studies dealing 
with phenomena such as radiation transport, atmospheric circulation, weather, and compo- 
sition should take aerosol effects into account (Egan et al^ 1977). 
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Stratospheric Observations 


The physical properties of the stratosphere can be probed using the Moon as a source 
as was done for six flights (Fig. 43-50) over the two-year period 1976-1977. Methane (3018 
cni~^), CO„ (3500-3800 cm~^), H„0 (3750-4000 cm”^) and other molecules can be identified 
from the spectra. The total column density of methane was foimd to be 5.8 x 10 molecules/ 

9 

cm“ above 14 km on 20 October 1977 (Smith and Hilgeman, 1978). This is equivalent to a 
niean mixing ratio of 1.9 ppmv, a value about 30% higher than other investigators have found for 
for the daytime stratosphere. Figures 43-47 were obtained with a beam size of 6.3 arc-min 
which has a limiting resolution of 10 cm~^. The resolution shown on Fig. 47 (R=1.8 cm”^) 
is not the actual resolution but what would be obtamed if the beam size were reduced. 

Figures 48-50 were measured wath a beam of 1 arc-min which corresponds to a limiting 
resolution of 2 cm A blocking filter with long wavelength cutoff of 3.15 M.m (3170 cm ^) 
was used on Fit. 20 (Fig. 48) preventing the strong CH^ bands centered at 3018 cm ^ 
from appearing. This cutoff filter was used for the M42 data of Fit. 19 (Figs. 37 and 38) 
and could not be replaced in time for Fit. 20. 

Examples of Reduction to Absolute Flux 

A discussion will now be given for converting the above sjxjctra which are shown in 
relative uiiits to absolute flux m Janskys (1 Jy=10 ‘'^W/m“/Hz) and true (labomtory) fre- 
quencies. The frequency can be corrected for by measuring lines of known frequency in the 
observed spectra, i.e. telluric II 2 O and ClI^ and plotting [v (true) - v (meas)] vs v (meas). 

This is done for the » Ori spectrum of Fit. 21 (Fig. 9) and the Moon spectrum of Fit. 

22 (File M00N22L1, Fig. 49) both at a resolution of 9.9 cm~^. The results ai’e shown m 
Fig. 51 where it is seen that the change between v (true) and v (meas) is about one-half of a 
resolution element (5 cm"^). The least square best fits from Fig. 51 ai*e now applied to the 
spectra of a Ori and the Moon to correct the frequencies. 

The spectrum of a Ori can now be divided by that of the Moon where the lunar sp)ectrum 

is assumed to have both a thermal component (303°K) and a reflected component (6000°K, 

4 -1 

albedo .01). The absolute flux of a Ori is normalized to 2.59 x 10 Jy at 4545 cm (K mag ~ 
-4.0) and is shown in Fig. 52. The strong depression of the continuuni from 4000-4400 cm”^ 
is caused by the first overtone brand of CO. Ideally, both the I'esponse of the telescope- 
interferometer system as well as the effect of the telluric absorption lines should divide out. 
hispeciicn of Fig. 52 shows residual CO„ bands in the region 4800-5200 cm” , as well as 
strong absorption in the 3550-3800 cm region which has been blocked out. The noise in the 
spectrum from 5200-5800 cm”^ is produced by the small value of the flux of both a Ori and 
the Moon in this region. 
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Figure 53 shows the results of a similar reduction fora Tau observed on 6 Feb 1976 
(Fit 7) with spectral resolution 13.3 cm"^. The region of the hot CO bands is shown with 
resolution 7.3 cm”^ (sine apodization). Line identification is facilitated with sine apodization, 
though the line shapes are distorted. The first overtone bands of and are 

shown. 

Mars is shown in absolute units in Fig. 54 for 11 Feb 1976 (Fit 9) with resolution 13.3 

cm"^. The CO„ bands at 4800-5100 cm~^, H^O bands at 3800-4000 cm~^ and water of 

-1 “ 

hydration bands 2900-3750 cm are shown. A line by line atmospheric modeling program 
has provided an analysis of these bands (Egan, et al , 1977, 1978). This modeling is necessary 
as the vibration-rotation lines of COg are saturated m both the telluric and Martian spectra, 
and therefore add m a nonlinear way. 

Figure 55 is a spectrum of Mira (o Cet) for 12 Feb 1976 (Fit. 10) with resolution 
13.3 cm~^ (solid lines). A theoretical model by Auman (1969) (dashed curve) with Te = 

2000°K, log g = -2 has been fit to the observed data and show's good agreement w'ith the hot 
steam bands centered at 3700 and 5300 cm but has too much flux in the region 3900- 
4500 cm This is to be expected as no CO band absorption w'as included in the model. 

SUMMARY 

During the period reported on, we built our optical system, developed a new and unique 
computer capability, and successfully took data on every flight scheduled. The data taken 
represent a scientific return an order of magnitude better than anticipated and should provide 
a fertile source for theoretical analysis for years to come. 

The excellent performance of the optical system served to validate our computer-aided 
optical design capability and demonstrated the predicted wavelength resolution limit of 
2.1 cm”'*’. In addition, the image quality for visual acquisition and guidance has been 
sufficient to enable us to observe stars dowii to V = 12 magnitude, an important consideration 
for some of the very red stars we have measured. 

Our new computer capability has proven its worth by allowmg us to make flight program 
decisions between each flight. The flexibility of its design became apparent when we had 
to reconfigui’c it to handle an order of magnitude greater data load than we had originally 
anticipated. 

Another important measure of performance is reliability, and this is especially true 
for linuted observing time programs. We were ready to take data for each scheduled flight, 
thus avoiding any wasted flight time. Our data Improved ^vlth each fllf^t in a sci'ies 
(sec Table 1) as we became aware of the potential of our equipment and kept "pushuig" it 
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further. Surprisingly little modification to our basic design was needed to achieve the 
order of magnitude improvement we eventually obtained. 

We have presented some of the better spectra obtained, illustrating spectra of 
standard stars, M stars, a nebula, planets and the moon. The spectra show an unprecedented 
capability to view at good spectral resolution through the hitherto opaque (near-ER) water 
bands. Water has been seen in the M stars, and in the surface of Mars. The full significance 
of these observations will take time to develop, but we have clearly demonstrated the 
capability to acquire data never seen before. 
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Fig. 1 Grumman FTS System Installed in NASA Lear Jet. 
Grumman Designed and Constructed Interferometer 
and Telescope System Shown as It Appeared After 
Flight Qualification Tests. Telescope Looks Through 
the Airplane Window (Not Shown) by Means of a 
Reflection Off the Airplane Motion Compensator. 
Entire Optical System Sits on a Vibration Isolated 
Platform 



Fig. 2 Schematic Optical System of Telescope and Interferometer 
as Installed in NA^A Lear Airborne Observatory 

ly 











Fig. 3 Optical Design Ray Trace Diagram Shows Path of Travel of Extreme and Center Rays of f/23 Beam as They 
are Reflected from a Field Mirror, a Collimating Mirror, the Interferometer, a Focusing Mirror, and a Folding 
Mirror to Reach the Detector Lens. Abscissa and Ordinate are Distances in Inches (Different Scales). Computer 
Aided Optical Design of This Type Allows Simulation of Proposed Systems Before they are Built 



Fig. 4 Completed System Ready for Testing. The 
Grumman Designed and Built Optical System 


Indicated in Figure 3 is Shown Here in a Com- 


pleted State Ready for Laboratory Testing 
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Fig. 5 Lunar Spectrum Taken from Ground at Grumman. Test Spectrum 
Shows Earth's Atmospheric Absorption Features Due to CO 2 end 
Water. Plot Shows Intensity of Lunar Signal (Uncorrected for Sys- 
tem Response) Versus Wavenumber in Inverse cm. Strong Water 
Absorption Makes Ground-based Astronomical Observations of 
Water Difficult 



Fig. 6a Spectrum of Alpha Ori on 5 Feb 1976 in 28004400 cm'^ Region at 13.3 cm'^ 
Resolution. 
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RELATIVE UNITS 


ALPHA ORI 



Fig. 7a Spectrum of Alpha Ori on 1 Sept 1976 in 2800-3600 cm'^ Region at 
1.8 cm'^ Resolution. 
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RELATIVE UNITS 


ALPHA ORI 



Fig. 7b Spectrum of Alpha Ori on 1 Sept 1976 in 3600-4400 cm'^ Region at 
1.8 cm'^ Resolution. 



Fig. 7c Spectrum of Alpha Ori on 1 Sept 1976 in 4400-5200 cm'^ Region at 


1.8 cm Resolution. 


23 








RELATIVE UNITS I RELATIVE UNITS 


ALPHA ORI 



Fig. 7d Spectrum of Alpha Ori on 1 Sept 1976 in 5200-6000 Region at 1.8 
cm'^ Resolution. 



Fig. 8a Spectrum of Alpha Ori on 1 Sept 1976 in 2800-3600 cm'^ Region 
at 2.1 cm'^ Resolution. 


2-t 





relative units 



Fig. 8b Spectrum of Alpha Ori on 1 Sept 1976 in 3600-4400 cm' ' Region at 
2.1 cm'^ Resolution. 


Fig. 8c Spectrum of Alpha Oii on 1 Sept 1976 in 4400-5200 cm' Region at 
2.1 cm'^ Resolution. 
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RELATIVE UNITS 


Fig. 8d Spectrum of Alpha Ori on 1 Sept 1976 in 5200-6000 cm'^ Region at 
2.1 cm'^ Resolution. 



Fig. 9a Spectrum of Alpha Ori on 8 Feb 1977 in 2800-4400 cm'^ Region at 
9.9 cm'^ Resolution. 
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RELATIVE UNITS 



Fig. 9b Spectrum of Alpha Ori on 8 Feb 1977 in 4400*6000 cm~^ Region at 
9.9 ccn*^ Resolution. 
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RELATIVE UNITS 



Fig. 10b Spectrum of Alpha Ori on 12 Oct 1977 in 4400-6000 cm'** Region 
at 9.9 cm'^ Resolution. 



Fig. 11a Spectrum of Alpha Ori on 13 Oct 1977 in 2800-4400 cm'^ Region 
at 9.9 cm'^ Resolution. 
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RELATIVE UNITS | | RELATIVE UNITS 



Fig. 11b Spectrum of Alpha Ori on 13 Oct 1977 in 4400-6000 cm'^ Region 
at 9.9 cm’^ Resolution. 



Fig. 12 Spectrum of Alpha Ori on 26 Oct 1977 in 2800-6000 cm'^ Region at 
37.7 cm'l Resolution. 
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Fig. 14a Spectrum of Alpha Tau on 6 Feb 1976 in 28004400 cm'^ Region 
at 13.3 cm'^ Reiolution. 
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RELATIVE UNITS 











Fig. 15b Spectrum of Alpha Boo on 27 Aug 1976 in 3600-4400 cm'^ Region 
at 1.8 cm'^ Resolution. 
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RELATIVE UNITS 



Fig. 16a Spectrum of Alpha Boo on 27 Aug 1976 in 2800-3600 cm'^ Region 
at 2.1 cm'^ Resolution. 
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Fig. 16b Spectrum of Alpha Boo on 27 Aug 1976 in 3600-4400 cm''* Region 
at 2.1 cm'^ Resolution. 



Fig. 16c Spectrum of Alpha Boo on 27 Aug 1976 in 4400-5200 cm*' Region 
at 2.1 cm'^ Resolution. 
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at 29.9 cm'^ Reiolution. 







Fig. 18 Spectrum of Mu Cep on 18 Oct 1977 in 2800-6000 cm*^ Region 
at 29.9 cm'^ Resolution. 



Fig. 19 Spect'"v. of Alpha Lyr on 9 Sept 1976 in 2800-6000 cm'^ Region 
at 56.0 cm'^ Resolution. 
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RELATIVE UNITS 




Fig. 21 Spectrum of Alpha CMa on 9 Feb 1977 in 2800-6000 cm'^ Region 
at 29.9 cm'^ Resolution. 
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RELATIVE UNITS 


ALPHA CM A 





Fig. 23 Spectrum of Alpha CMa on 18 Oct 1977 in 2800-6000 cm'^ Region 
at 37.7 cm'^ Resolution. 








Fig. 24 Spectrum of RX Boo on 6 Feb 1976 in 2800-6000 cm*^ Region 
at 29.9 cm’^ Resolution. 



Fig. 25 Spectrum of RX Boo on 27 Aug 1976 in 2800-6000 cm ^ Region 
at 29.9 cm'^ Resolution. 
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Figure 26a. Spectrum of Mira on 12 Feb 1976 in 2800-4400 cm'1 Region 
at 13.3 cm'1 Resolution. 



Fig. 26b Spectrum of Mira on 12 Feb 1976 in 4400-6000 cm'^ Region 


at 13.3 cm'^ Resolution. 
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Fig. 27a Spectrum of Mira on 10 Feb 1977 in 2800-4400 cm‘^ Region 
at 9.9 cm’^ Resolution. 
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Fig. 29 Spectrum of Mira on 14 Oct 1977 in 2800-6000 cm’^ Region 


-1 


at 29.9 cm Resolution. 



Fig. 30a Spectrum of R Leo on 17 Feb 1977 in 2800-4400 cm'^ Region 
at 9.9 cm'^ Resolution. 


13 











Fig. 30b Spectrum of R Leo on 17 Feb 1977 in 4400-6000 cm"^ Region 
at 9.9 cm'^ Resolution. 



Fig. 31a Spectrum of R Leo on 18 Feb 1977 in 2800-4400 cm ^ Region 
at9.9cm'^ Resolution. 
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Fig. 31b Spectrum of R Leo on 18 Feb 1977 in 4400-6000 cm’^ Region 
at 9.9 cm*^ Resolution. 



Fig. 32a Spectrum of R Leo on 18 Feb 1977 in 2800-4400 cm'^ Region 
at 9.9 cm'^ Resolution. 
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RELATIVE UNITS 



Fig. 32b Spectrum of R Leo on 18 Feb 1977 in 4400-6000 cm ^ Region 
at 9.9 cm'^ Resolution. 



Fig. 33a Spectrum of R Leo on 18 Feb 1977 in 2800-4400 cm'^ Region 
at 13.3 cm'^ Resolution. 
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Fig. 33b Spectrum of R Leo on 18 Feb 1977 in 4400-6000 cm'^ Region 

.*1 

at 13.3 cm Resolution. 



Fig. 34 Spectrum of R Leo on 18 Oct 1977 in 2800-6000 cm'^ Region 
at 37.7 cm'^ Resolution. 
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Fig. 35 Spectrum of U Her on 19 Oct 1977 in 2800-6000 cm’^ Region 
at 37.5 cm'^ Resolution. 



Fig. 36 Spectrum of R Aql on 20 Oct 1977 in 2800-6000 cm'^ Region 
at 37.3 cm'^ Resolution. 
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Fig. 37 Spectrum of M 42 on 9 Sept 1976 in 2800-6000 cm''^ Region 

at 56.0 cm'‘l Resolution. The Hydrogen Emission Lines Paschen a, 
Brackett 0, and Brackett y are Shown. 



Fig. 38 Spectrum of M 42 on 9 Sept 1376 in 2800-6000 rm*^ Region at 
56.0 cm’1 Resolution in Janskys. The Hydrogen Emission Lines, 
Patches a, Brackett Q and Brackett y are Shown. 
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Fig. 39 Spectrum of M 42 on 27 Oct 1977 in 2800-6000 cm'^ Region at 
56.0 cm'‘l Resolution. The Hydrogen Emission Lines Paschen <a 
and Brackett y are Shown. 



Fig. 40a Spectrum of \/enus on 17 Feb 1977 in 28004400 cm’^ Region 
at 9.9 cm'^ Resolution. 
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RELATIVE UNITS 



Fig. 40b Spectrum of Venus on 17 Feb 1977 in 4400-6000 cm‘^ Region 
at 9.9 cm‘1 Resolution. 



Fig. 41a Spectrum of Mars on 11 f eb 1976 in 2800-4400 cm'^ Region 
at 13.3 cm*^ Resolution. 
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Fig. 42 Spectrum ot Mars on 28 Oct 1977 in 2800-6000 cm''l Region 
at 29.9 cm'^ Resolution. 
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Fig. 43a Spectrum of Moon on 10 Feb 1976 in 2800-4400 cm’^ Region 
at 13.3 cm'1 Resolution. 



Fig. 43b Spectrum of Moon on 10 Feb 1976 in 4400-6000 cm'^ Region 
at 13.3 cm*'l Resolution. 
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Fig. 45b Spectrum of Moon on 11 Feb 1976 in 4400-6000 Region 
at 13.3 cm'^ Resolution. 









-1 


Fig. 46a Spectrum of Moon on 11 Feb 1976 in 2800-4400 cm Region 


at 13.3 cm' Resolution. 











Fig. 47a,b Spectrum of Moon on 3 Sept 1976 in 2800-3600 cm-1 region 
at 1.8 cm-1 resolution. Upper curve (a) is L aperture and lower 
curve (b) is R aperture. 



Fig. 47c, d Spectrum of Moon on 3 Sept 1976 in 3600-4400 cm-1 region 
at 1.8 cm-1 resolution. Upper curve (c) is L aperture and lower 
curve (d) is R aperture. 
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Fig. 47e,f Sp«ctrum of Moon on 3 Sept 1976 in 5200-6000 cm' region 
at 1.8 cm-1 resolution. Upper curve (e) is L aperture and low- 
er curve (f ) is R aperture. 



Fig. 47g,h Spectrum of Moon on 3 Sept 1976 in 4400-6000 cm-1 region 
at 1.8 cm-1 resolution. Upper curve (g) is L aperture and low- 
er curve (h) is R aperture. 

58 









Fig. 48a, b Spectrum of Moon on 10 Sept 1976 in 2800-3600 cm-1 region 
at 2.1 cm-1 resolution. Upper curve (a) is L aperture and low- 
er curve (b) is R aperture. 



Fig. 48c,d Spectrum of Moon on 10 Sept 1976 in 3600-4400 cm-1 region 
at 2.1 cm-1 resolution. Upper curve (c) is L aperture and low- 
er curve (d) is R aperture. 
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Fig. 48e,f Spectrum of Moon on 10 Sept 1976 in 4400-5200 cm-1 region 
at 2.1 cm-1 resolution. Upper curve (e) is L aperture and lower 
curve (f) is R aperture. 



Fig. 48g,h Spectrum of Moon on 10 Sept 1976 in 5200-6000 cm-1 region 
at 2.1 cm-1 resolution. Upper curve (g) is L aperture and lower 
curve (h) is R aperture. 
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Fig. 49a, b Spectrum of Moon on 9 Feb 1977 in 2800-3600 cm-1 region 
at 2.1 cm-1 resolution. Upper (a) is L aperture and low- 
curve (b) is R aperture. 



Fig. 49c,d Spectrum of Moon on 9 Feb 1977 in 3600-4400 cm-1 region 
at 2.1 cm-1 resolution. Upper curve (c) it L aperture and low- 
er curve (d) it R aperture. 
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Fig. 49e,f Spectrum of Moon on 9 Feb 1977 n 4400-5200 cm-1 region 
at 2.1 cm-1 reiolution. Upper cuive (e) is L aperture and low- 


er curve (f| is R aperture. 



Fig. 49g,h Spectrum of Moon on 9 Feb 1977 in 5200-6000 cm-1 region 
at 2.1 cm-1 resolution. Upper curve |g) is L aperture and low- 
er curve (h) is R sparture. 


()2 







Fig. 50a i'ntctrum of Moon on 20 Oct 1977 in 2800-3600 cm-1 region 
a* VU cm-1 resolution 



Fig. 50b Spectrum of Moon on 20 Oct 1977 in 3600-4400 cm-1 region 
at 1 .8 cm-1 resolution 
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Fig. 50c Spectrum of Moon on 20 Oct 1977 in 4400-5200 cni-1 region 
at 1.8 cm-1 resolution. 



Fig. 50d Spectrum of Moon on 20 Oct 1977 in 5200-6000 cm-1 region 
at 1.8 cm-1 resolution. 
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Fig. 51 Frequency Calibration of Alpha Ori and Moon. A Least Squares 
Straight Line Has Been Fit to the Data Which is Used to Correct 
the Measured Frequencies in the Spectra r (Meas) to Laboratory 
Frequencies)' (True). 



Fig. 5? Spectrum of Alpha Ori on 8 Feb 1977 in 2800-6000 cm'^ Region 

at 9.9 cm’^ Resolution in Janskys (10'^® W/m^/H/). Tlic CO Lines 
at 4000-4360 cm’^ Are Shown. 
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Fig. 53 CO Bands of a Tau With Spectral Resolution 7.3 cm- ' (Sine Apodization). The First 
Overtone Bands of ^®0 And ^®0 Are Shown 


MAR5 11 TEB 137B 


WAVENUMBER <CM-1> 


Fig. 54 Spectrum Of Mars, 11 Feb 1976 With Spectral Resolution 1 3.3cm~V CO 2 Bands At 4800-5100 
H 2 O Bands At 3800-4000 cm'^. And Water Of Hydration Bands 2900-3750 cm-1 Are Shown 







MIR^ 12 "TB igVB 



Fig. 55 Mira (o Cet), 12 Feb. 1976 With Spectral Resolution 13.3 cm'^ (Solid Line). A Theoretical Model 
Of Auman, 1969 (Dashed Curve) Has Been Fit To Observed Data And Shows Good Agreement With Hot Water 
Bands Centered at 3700 and 5300 cm'\ But The Model Predicts Too Much Flux In Region Of 3900-4500 cm'V 
This is To Be expected as No CO Band Absorption Was Included In Model 
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APPENDIX 


BASIC programs, PHASE, and HFFT used in the data reduction procedure are 
listed. Some comments are included. 
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^\sE 

10 NI IS NO. OF HUINfS I M 1 N T K'<F t ^DG ^ A (MAX = 7100 ) 

20 -ItM N IS NO. OF POI-'^S JSEJ 10 GKNt.lAFF PHASES 
30 ’EM N 9 IS f'OSIflON OF SYMlETin POINT OF I NT E ^FE'^OGPAh 
40 HEM K 9 IS NO, OF I’MAbE CO^hl, OESMEO 
50 '^EM N 7 IS length of CONVO.VING FCTN. 

GO JIM XI 7100 ltl( 26 ]iS(‘'.?l.«(o 2 )«O(S 0 ) 

70 JIM Aslai.Wblft) 

DO 'l :50 
90 N7-51 
100 X 9=3 

110 files *,* 

120 irU’Uf "INI’OT INILKI- El<OGPA^ FILE NAME "lAS 
130 PNI jT 

140 ASSIGN A". , 1 »S 0 

150 IF so-u Then ido 

160 POINT "SOMETHING ..KONG Kill FILE NAME " 

170 STOP 

180 INPUT "InT'OT no, of PUInTS IN I N f E <F£ LOGJAM ",Nl 
no PKINT 

POO INPUl "OUTPUT PHASE COKHECIEJ HLE N.\ IF (NO N-JEJ (0 CKEAT*) " , a (. 
210 I’KINI 

220 ASSIGN • 3 'f, 2 ,K 0 
230 IF KU =0 THEN 270 
■VtO IF K 0<>3 T'IcN 160 
2,0 ;::.'i;a(F hi ,.I".,SS. 2 S 6 
. '.II IF 21 0 f.:L N 220 

2/0 l.i.’Ul "iNl’iJf IN .iNCMi.UJ:, JAM SYMMETHi' HJINT ",N9 
2'l,i i'HInT 

290 01 FNO »1 THEN 310 
100 mat kEAJ ,',‘ 1 ;X 


110 <, a 9*1 

32 o -Ni .,N/*i 
130 N 2 =(.N/ 2)*1 
340 N 3 =N /2 
350 N 5=9 

360 FOW XI=l TO K 9 
370 FOK 1=1 TO N 

380 S( I ]=( 6 . 2 DJ 1 D/N)«(I-L) 

390 NEXT I 
400 N 8 -:N 9-1 

410 FOH 1=1 TO NJ 

420 UimxlND*!) 

430 UlN»l-( l=X(N‘m 1 

440 NEXT I 
450 PPInT 01 I J 
460 PHInT 
470 FOH J=NS TO N 2 
480 H(J]=I(JJ =0 

490 KJH <=1 TO N 

500 d=(J-l)*S(K) 

510 H(J] = H( jJ.OIKlttCOSO) 

520 II J)=I(J)- 01 K)»SIN( 3 ) 

530 NEXT K 

540 EXT J 

550 F=7901 . 14*2/N 

560 FOH L=N 5 *l TO N 2-1 

5 70 0=.54*HIL J^.^J^lHIL-l J*H(L*l 1) 

580 t=.54*IlU*.23*(IIL-l ]*IlLM J) 

590 P= (ATN(E/0) ) *57.2953 
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'lOO IF D<0 iHtN f)20 

(Sin GOTO 630 

620 P=P»lfiO 

630 IF P>1 J THEN 650 

660 GOTO 660 

650 P=P-360 

660 U(LJ=P 

670 NEAT L 

680 Q1N21=0 

500 ('PINT " INPUT ;>MA5E5" 

700 PPINT 
710 05=0 

720 FOW I=N5*1 TO M2 
no G=Fo(I-l) 

760 point ItG.Om 

750 U5=05*AUS(UI IJ) 

760 0m-0[ I J/57.2058 

770 S(2“I-1 l=C0b(0n J) 

790 bI2*n = SlN(tl[ I 1) 

790 NEXT I 
000 PI71NT 

HO print "AVERAGE ABSOLUTE JIKFtRENCE - " i Tj/ ( N2-N5) 

820 i6<InT 

no A1=6.28319/N 

960 FOR J=1 TO n7 

950 Ml- J- ( (N7 -1 ) /2) -1 

860 M2: Ip* J. l6lS9/( (N7 1 ) /2) 

870 Sl:S2-:0 

880 FOR I-N5>1 fO J2 

890 A2 :A1<* ( I-l ) «(A1 

000 SI -S1*SI2«I-1 ]-'COS(A2) 

010 S2 = S2»S(2<*1 )*S' J( A2) 

020 NEXT I 

930 Ul=5l-S2 

960 RI J1 = Q1'*(.56».';5-*CUS(-12) ) 

950 RlJl=RtJ)/N2 

960 NEXT J 
970 print 

980 PRINT " C 'V.lLVrO Ij TL (KE8)'i < Ai" 

900 print 
OOO FUR J=1 TO :19 

,010 A=0 

1020 FOR 1=1 TO N7 

1030 A = A*Kl I 1<*XI I ♦J-l J 

1060 NEAT I 

1050 XtJlsA 

1060 IF <1=K9 THEN 1080 

1070 GOTO 1090 

1090 PRINT f«2lXlJ] 

1090 NEXT J 
1100 PRINT 

1110 print " (JIFFtRFiNCE (lE^T - "IGID" 

1120 print 

1130 8H=N9- ( (N7-1 ) /2) 

1160 FOR J=1 TO 11 

1150 K=J-1 

1160 L = X[88-n 

1170 R=X(88*A) 

1190 C=L-R 

1 190 U( J}=A8S(C) 
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l.’OO HHIMT K.LtW.C 

l-’lO ^JEAf J 
J = 0 

I ? )0 FOM J=1 rO 11 

U’'tO U = 0*(J(JJ 

12Sn MEAT J 
U’rjn 0-0/ 10 
l?/n t’KlNT 

I-’KImT "AVEKAfjE AilSULlHE JIFFEUEMCE = "JO 
I’KIMT 
1.100 I’KlMT 
1310 N9-M8 

13^n M9-’10’M7*1 

rn.) MEAT K1 

) J-Hl "i’OSIlION UF 'iY-l-lElMY POINT = "(Ml 

115D ,’<INT 

llr,(l PPIMT "CPU bECUMDS = 'MCPJ(O) 
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HFFT 

m FFT does cosine IRANSfORM ONLY 
?n REM USING HAMMING Ai’UOUAriON FCfN 

30 MEM N7 IS roiAL NO. OF DATA RTS FROM InTERF SYMM PT 
40 REM N9 IS InTEkFERUGRAM SYMMEFRY POINT 

50 REM G determines NO, OF PTS TO HE TRANSFORMED (N = 2-nG) 

60 REM W1 AND Wl ARE UAVENUM3ER LIMITS DESIRED (MIN.MAA = 2800, SOOO) 
70 DIM M4100],Y(4100J9A[4,40I 
80 FILES <»#» 

90 DIM AM81,lltl8J 

loo INPUT "INPUT PHASE CORR, I N TCRF E.'OSRA M FILE NAME' ",AS 
110 PRINT 

120 ASSIGN A«;,1,S0 
130 IF S0=0 then ISO 
140 IF S0<>3 THEN 100 

150 INPUT "INPUT TOTAL NO. OF 3ATA PTS KRDm INTERF SYNM PT ",n7 
160 PRINT 

170 INPUT "Input interf symmetry pt ",no 
180 PRINT 

190 INPUT "OUTPUT SPECTRUM F I lE NAME (NO NFED TD CREATE) ",Hk 
200 PRINT 
'10 ASSIGN •J^,2»R0 
■20 IF t<0-0 THEN 2 70 
230 IF H0<>3 THEN 190 
240 C>iEATE Ml ,HS, 130»2S6 
250 IF Rl-0 THEN 210 
260 srup 

2 70 rr’uT "nror value of g (no. of pis ro ie T'T-.ns" = 2aG) ",g 
200 print 

290 N-2»0G 

300 RG=l.Hli«7901.14/N7 

310 INPUT "INPUT .■■■AVENU-IDER LIMITS (./l,i,'2) ",wl,./2 

320 PRINT 

330 F=7'901.14/N 

340 J3=INT((Wl/F)*l) 

350 J4-INT ( (W2/E) ♦!) 

360 mat X-2ER 

3 70 MAT Yi^XER 

.300 FUR 1 = 1 TO N'J-l 

390 read ;<i«a 

400 NEAT I 

410 FUR 1=1 TU N7/2 

420 read «1 lAl I*2)»Y( I*2J 

430 NEAT I 

440 RF:aD rtHXIN/2»3] 

450 Y(N/2*31=0 

460 S1=S2=0 

470 FOR 1=1 TO N/2*l 

480 S1=S1*A(I*21 

490 S2-:S2*Y( I«2) 

Soo NEAT I 

510 bl-2'*51“AUJ-AlN/2*3) 

520 S2 = 2<»S2 
530 AI 1 l=xl 31 
540 Y(1)=0 

550 FUR 1=2 TO N/2*l 
560 A(11=A(I«2) 

570 Y( 1 ) = (Y( I*2J-YI 1*1 J)/2 
580 NEAT I 
590 Y(N/2MJ=0 
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ORfGlNAL PAQi |S 
OF POOR QUALfTt 



60n T = 2»3, 1'*l59/N 
610 M2=N/2*2 
520 'J3 = N/2*J 
530 fOW J=1 TO N/6*l 
5(»0 J1 = J 
550 FOH L=1 TO 2 
'-50 n = To(j-l) 

b70 K-:N/2*2-J 

690 l.l='J/2*l*L 

590 S-SI»4(m 

700 C=;C0S(T1) 

710 X1=X(J]*X(K] 

720 X2=X(JJ-X(KJ 

730 n=riJ)*Y(X) 

740 Y?=YtJl-Y[X) 

750 X(L1 )=<Xl«X2*S-Yl»C)/2 

750 Y[L1] = <Y2»Y1‘*S*X2"C)/2 

770 J=K 

790 next L 

790 J=J1 

900 i<=N2-J 

910 X(J]=X(N21 

920 Y[J)=Y(N2) 

930 X(M = X(N3) 

890 Y(M-Y(N31 

150 NEAT J 
I ') 0 5 " G •• 1 
ii70 let 

990 let 9=U«ATN(1)/N 

990 EOW 1=1 TO N 

900 LEI xm=XII)/N 

910 LET Ym = Ym/NI 

920 NEXT I 

)30 FOh< L = 0 TO G-1 

9!f0 LET Gl=2<** (G-L-1 ) 

950 LET M = 0 

950 FOK 1=1 TO 2*®L 

970 let Kl-INTIM/Gl) 

980 1690 

990 LEI Y1-C0S(2'>K2> 

1 000 let Y2 = -SIN(^^'>K2) 

1010 FOrt J=1 TO G1 

1020 Ml=MfGl*l 

1 030 r-l2 = M*l 

lOOO Y3 = X( 81 J Yl-YCMl J<*Y2 

1J5.1 Y4 = XIM1 J*Y2*Y(M1 1»Y1 

1050 X[M1 ) = X[M2 1-Y3 

1070 YlMl J=Y182 J-Y4 

1090 X182]=X(M2 J*Y3 

1090 Y[M21=Y(82'*Y4 

1100 LET M=8*l 

1110 NEXT J 

1120 let M=M*Gl 

1130 NEXT I 

1140 MEXT L 

1150 FO^ 1=0 TO N-1 

1150 LET M = I 

1170 GObUrf 1690 

1190 IF K2>rl IHEN 1270 

1190 I1=I*1 
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i?on K4-K2*i 
12ln (<3 = x(m 
l?2n x(in=x(K<*] 

1?30 X(K<»1:K3 

1240 K3=r(ii) 

I25n rtIlJ=Y(K4) 

1260 Y(K6J=K3 
1270 MEXT I 

IP 10 X( 1 l=(Sl*S2)/(0 >N) 

*290 ri = 3. 16159/(2'*N) 

1300 'U=N*1 

1310 'J2-N*2 

1320 fOK J=l TO N/2 

1 130 J2 = J 

13'm) ^ OK L-1 TO 2 

1350 J1^J»2-L 

1360 K=N1-J 

1 370 xi = 'J2-Jl 

1330 Lls^J^L 

1390 A1=XC J1 J*X[Kl J 

1000 A:Y(JJ*Y(Kl 

r»10 t)l-(Xl J1 1-X(K1 l)/SlM(n*(2^Jl-2) ) 

1620 u=(yuj-y(6])/sin( ri^ca^j-i) ) 

1630 A(L1 l-(Al-Ul>/6 

1660 YlLl ]=(A -d)/4 

1 650 J :K 

l':60 JEAT L 
1 6 70 JvJ2 
16 10 6-li*l-J 
1 690 XI J*1 )=X[ 'Jl J 
1 500 Y[J] = Yl.Nl] 

1510 AllX]--X(ri2J 
1520 Y(r. 1-YIM21 
1530 NEXT J 
1560 31=11 
1550 II=INT(J3/2)*1 
1560 12=INT (J6/2) ♦! 

1570 L1=I2-I 1»1*2»‘U 

1530 fOK Irl TU LI 

1590 X[6101-n=xn2M-I*31 ) 

1600 NEAT I 

1610 POK I=I1-M1 TO 12*31 
1620 X = 2» ( I-l 1*M1 ) *1 

1630 XIM = X[6100-I2*I- II J 

1660 XlK*l)rYm 

1650 NEXT I 
1660 X=1 

1670 J5=2* ( I 1-Ml ) -1 
1680 30r0 1 770 
1690 let ^2=0 
1 700 ^'OK 6sl TO 5 
1710 lET »X3=INr (Kl/2) 

1720 LEI K2=2« (K2-K31 *X1 

1710 LET K1=k3 

1760 NEXT K 
1750 KETUKM 

1710 FOLLOWING InTE ISOLATES ANO SMOOTHS FOW COSINE T46NSF0K3 OnlY 

17/0 N3=20 

1730 3s2«N8 

1790 POK K=1 TO 4 
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I'lon 

run 

IT30 

in<*o 

lfl50 

1H60 

run 

I'i'lO 

1!V)0 

r>on 

1910 

io?n 

1910 

rr»n 

19 S 0 
1950 
1970 
1910 
1990 
flOoo 
/?0 1 o 
i'Oi'O 
^010 
/^O'tO 
-’() iO 
/050 


Fun I=l ro M 

5= 1. i'*rrr» (N , i ♦...’■<) 

U= (3. l'♦l'j9/^^) 

AlKf I r-. (SIN (15) / AM ( ,": , - .5 5/ ( 1 -C) ) 

Nf.XT I 
NtXI < 

J3-J3-J5M 
J5=J5-J3*1 
Fun j:J3 TO J5 

n 1 1 - ,S5'^x( 1 ♦ ./ j > (.<( j 1 1 'A( jM n 

Y(«;i = Y(3)-Y( l-Y(Sl-0 

FOK K=1 ro 5 
Fon I-l 10 M 

Y( * 1 )-YUM ] t.Yf M ( • ru MK, I) 
nf.xt I 
NtAf < 

FOK 1=1 TO 5 

print o^JYII) 
next I 
NEXT J 

I’KINI "CUC. '..'AVENU'iilEW RINSE " ; I" •» ( J i ^ V ) « " I U” ; F ■> ( JA ♦ J5-6/5 ) 

print 

print "riAVENUMuER INrEI<VAL 'UF/5 
print 

print "RESOlU ion = "5R6»"CN-1'' 

■'RIN r 

.■■NINI "CPU NFcO'lUS : 'UCPI(n) 


/ 
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